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We derive the generalized Ginzburg-Landau free energy functional for conventional and unconven- 
tional singlet superconductors in the presence of paramagnetic, orbital and impurity effects. Within 
the mean field theory, we determine the criterion for appearence of the non uniform (Fulde-Ferrell- 
Larkin-Ovchinnikov) superconducting state, with vortex lattice structure and additional modulation 
along the magnetic field. We also discuss the possible change of the order of transition from normal 
to superconducting state. We find that the superconducting phase diagram is very sensitive to 
geometrical effects such as the nature of the order parameter and the shape of the Fermi surface. In 
particular, we obtain the qualitative phase diagrams for three-dimensional isotropic s-wave supercon- 
ductors and in quasi two-dimensional d-wave superconductors under magnetic field perpendicular 
to the conducting layers. 

In addition, we determine the criterion for instability toward non uniform superconducting state 
in s-wave superconductors in the dirty limit. 
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I. INTRODUCTION 

The relative importance of orbital and paramagnetic 
effects in suppression of superconductivity is determined 
by the ratio of the orbital upper critical fields if C 20 ~ 
$o/27r£oa£o& and the paramagnetic limiting field^ H p = 
w 0.7lAo//i, called the Maki parameter— olm — 
V2H C 2o/H p . Here 4>o is the flux quantum, £oa and £ofc 
are the superconducting coherence lengths in two mutu- 
ally perpendicular and perpendicular to magnetic field 
directions a and &, A is the superconducting gap at zero 
temperature and fi = is the electron magnetic 

moment. Usually, the Maki parameter is of the order 
of the ratio of critical temperature to the Fermi energy 
olm ~ T c /sf- That demonstrates the negligibly small 
influence of paramagnetic effects on superconductivity. 
However, in the case of small Fermi velocity (that hap- 
pens in materials with heavy electronic effective mass) 
or in the layered metals under magnetic field parallel to 
the layers, the value of the Maki parameter can be even 
larger than unity. 

The consideration of the magnetic field acting only 
on the electron spins, corresponding to the limiting case 
of infinitely large Maki parameter, leads to some pecu- 
liar effects. First, the phase transition from the normal 
metal to the superconducting state which is of the sec- 
ond order in low field - high temperature region changes 
to the first orde»Si£ at fields above H* « 1.06T c //x 
and temperatures below T* sa 0.56T C . Starting at 
this critical point, the line of the first order transi- 
tion is finished at zero temperature and at the mag- 
netic field equal to the Chandrasekhar-Clogston limit- 
ing field Hp. However, as was shown by Fulde and 
Ferrelfi and Larkin and Ovchinniko\iiS, even at larger 
field i/FFLO ~ 0.755A /yU « l.07H p , the normal state 
is unstable with respect to the second order type tran- 
sition to the inhomogeneous cosine like gap modulated 



superconducting state (FFLO state) with wave vector 
q c w 2.38/i.ffFFLo/i'F- The recent calculations^! at zero 
temperature have demonstrated that more complicated 
crystal structures are more favorable than the simple 
plane wave. A first order type transition to the face- 
centered cube superconducting state was predicted to 
occur at the field larger than H^flo- This conclusion 
is in correspondence with the finite temperature inves- 
tigations performed in vicinity of critical point showing 
the appearence of FFLO superconducting state below the 
critical temperatur e) 12 : 13 . 

These results are changed a lot due to effects of orbital 
depairing and impurities. 

The role of the orbital effects was studied first at T = 
by Gruenberg and Gunthes^ who have demonstrated 
that the FFLO state appears in pure metal (assuming 
that it is formed by means the second order transition) 
if the Maki parameter is larger than 1.8. 

The influence of impurities in absence of the orbital ef- 
fect was investigated by Aslamazovii^. He found that im- 
purities do not kill the FFLO state but decrease the field 
-Hfflo °f absolute instability of the normal state for the 
FFLO formation such that, in the dirty limit (tT c <C 1), 
at zero temperature, -/Jfflo is lower than the field of 
the first order transition to homogeneous superconduct- 
ing state, Hp. Physically, it does not yet abolish a pos- 
sibility of existence of inhomogeneous superconducting 
state because the actual phase transition from the nor- 
mal state could be of the first order transition to FFLO 
state at some field H > H p . 

The investigation of orbital effects near the critical 
point was performed for isotropic three-dimensional pure 
metals by Houzet and Buzdini^. It was found that, unlike 
the conclusions obtained in the absence of orbital effect, 
for finite but large enough Maki parameter, the FFLO 
modulated state arises from the normal state starting 
from some temperature higher than the critical temper- 
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aturc. 

All the studies cited above concerned the case of 
isotropic s-wave superconductivity. The theoretical in- 
terest to the FFLO state in superconductors with d- 
painng was stimulated by the experimental 

identification of the pairing state in several of the high-T c 
cuprate superconductors and heavy fermionic materials. 

The recently discovered heavy fermionic tetragonal 
compound CeCoIn^ was established as d x 2_ y 2 super- 
conductor similar to high-T c cupratesS 1 * 2 ^. In this com- 
pound, the phase transition to the superconducting state 
becomes of the first order at low temperature - high field 
region and possible formation of FFLO at lower tem- 
peratures was reported for the magnetic field directed 
parallelSiMiS^ as well as perpendicular— to the basal 
plane. 

The first theoretical investigation of the phase diagram 
in the tetragonal, doped-by-impurities superconductor 
with d-pairing under the field parallel to c-axis was done 
in the Ref— . It was found that, in the absence of orbital 
effect, the change of the type of transition from the sec- 
ond to the first order occurs at some temperature which 
is lower than the temperature of appearence of FFLO 
state. 

The orbital effect in the same type of superconduc- 
tor with quasi two-dimensional spectrum was taken into 
account by Ikeda and Adachi21 and the different phase 
diagram topology was established. That is, in contrast 
with clean s-wave isotropic superconductor, the FFLO 
state arises from the normal state starting from some 
temperature lower than the critical temperature. This 
result was ascribed by the authors of Ref— to the non- 
pertubative treatment of the orbital effects incorporated 
there. 

It seems, however, that, in absence of analytical cal- 
culations, it is dificult to recognize an inequivocal reason 
for this discrepancy. The main goal of the present article 
is to make clear the influence of paramagnetic, orbital 
and impurity effects on the phase transition of a nor- 
mal metal to superconducting state, including the FFLO 
state formation and type of the phase transition. With 
this purpose, we shall derive Ginzburg-Landau functional 
for the conventional and unconventional superconduct- 
ing state with singlet pairing in the metal with arbitrary 
point symmetry and with arbitrary amount of point like 
(s-wave scattering) impurities. Then, for the cases of 
isotropic metal with s-pairing and tetragonal supercon- 
ductor with d-pairing under magnetic field parallel to c- 
axis, the simple analytic criteria of appearence of FFLO 
state and the type of normal-superconductor phase tran- 
sition shall be established. In particular, we shall demon- 
strate which temperature of FFLO appearence or the 
critical point temperature is higher. 

The structure of the article is as follows. We begin with 
the general expressions of Ginzburg-Landau functional 
for the superconducting state (in metal with the arbi- 
trary concentration of impurities) transforming accord- 
ing to identity and non-identity representations of the 



crystal point group symmetry. The corresponding deriva- 
tion of this functional from microscopic theory valid at 
finite temperature in vicinity of critical point is found 
in the Appendices. Then, for the cases of s-pairing and 
d-pairing, the criteria of FFLO state existence and the 
first order type transition and their competition shall be 
formulated. In addition to these finite temperature cal- 
culations, the critical field of dirty normal metal insta- 
bilty to FFLO state formation in presence of the orbital 
effect (generalization of the papers by Gruenberg and 
Gunther— and by Aslamazo\ii^) at zero temperature is 
found. 



II. FREE ENERGY NEAR CRITICAL POINT 

The Ginzburg-Landau functional consists of the sum 
of the leading terms in the expansion of the supercon- 
ducting free energy in the order parameter A and its 
gradients. 

In the purely orbital limit, it contains terms propor- 
tional to |A| 2 , |A| 4 and |VA| 2 , with coefficients de- 
pending on the temperature and impurity concentration. 
Strictly speaking, it is only valid near the critical tem- 
perature T c of the second order transition from normal 
to superconducting state, when the coefficient a in front 
of |A| 2 vanishes. Close to this point, the amplitude of 
the gap is indeed small and the magnetic length which 
determines the characteristic scale of the variation of A 
coincides with the thermal correlation length which di- 
verges at T c . This allows to retain only the first term in 
gradient expansion. 

In the purely paramagnetic limit, the coefficients in the 
functional also depend on the magnetic field. The equa- 
tion a(T,H) = then defines the transition line Hq(T) 
from normal state to uniform superconducting state in 
the temperature-magnetic field phase diagram. Along 
this line, the coefficient 7 in front of |VA| 2 happens to 
change its sign. This signals the instability toward the 
modulated FFLO superconducting state. In order to es- 
tablish the modulation wavelength in FFLO state, higher 
order terms in the gradient expansion should also be in- 
cluded in the functional. One can restrict the free energy 
expansion to the term | V 2 A| 2 only in the vicinity of the 
triple point, where the FFLO instability occurs and the 
typical FFLO modulation wavelength diverges. The co- 
efficient j3 in front of |A| 4 also may change its sign. This 
signals a critical point when the type of the transition 
into uniform superconducting state changes from second 
order to first order as the temperature is lowered. In such 
case, the type of the transition into FFLO state will be 
determined by the sign of the fourth order terms in A 
and of higher order in the gradient expansion. Again, 
close to the critical point, one can consider the terms of 
the order |A| 2 |VA| 2 onlyi£. The peculiarity of the mi- 
croscopic theory is that, in the pure limit, 7 and (3 change 
sign at the same place, with coinciding triple and critical 
point: the tricritical point (T*,H*), with H* = H (T*). 
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In the presence of impurities, however, the triple point 
and the critical point do not coincide any longer. For 
s-wave superconductors, the triple point occurs at lower 
temperature than the critical point, while for d-wave su- 
perconductors, the opposite situation takes placed. 

The effect of the orbital field on the interplay between 
transition into conventional superconducting state (with 
vortex lattice in such case) and FFLO state (with FFLO 
modulation in the direction parallel to the vortex axes) 
was considered within this frame in s-wave superconduc- 
tors and in pure limit onlji— As it is important that the 
magnetic length remains large compared to the supercon- 
ducting coherence length, a Ginzburg-Landau expansion 
is only possible when the paramagnetic effect is much 
larger than the orbital effect (large Maki parameter) . In 
the pure s-wave superconductor, it was shown that the 
triple point was moved to higher temperatures— There- 
fore, impurities and orbital effect act in opposite direc- 
tions in the s-wave case. 

The goal of the two next Sections is to provide a frame 
to discuss the nature of the transition from normal to 
superconducting state in the presence of impurities and 
a small orbital effect in superconductors with arbitrary 
(even and one-component) order parameter. 

The free energy Ginzburg-Landau functional up to 
the fourth order terms of the order parameter and the 
fourth order terms in gradients for the isotropic s-pairing 
superconductors doped by impurities has been derived 
first in the paper—. With a purpose of investigation 
of FFLO state in clean s-wave superconductors a simi- 
lar result based on calculation using Eilenberger— and 
Larkin and Ovchinnikov 30 formalism has been derived 
also in the paper—. The derivation for the dirty d-wave 
tetragonal superconductor based on the direct calcula- 
tion of the vertex parts renormalization first introduced 
by Gor'ko\i2i has been accomplished in the paper—, then 
including all orders in gradients in the paper—. Our 
derivation made for reliability by both Gor'ko\s2i and 
Eilcnbcrgei— i2£ methods (see Appendices A and B) is 
related to the case of doped-by-impuritics superconduct- 
ing metal of arbitrary crystaline symmetry with an order 
parameter 



A(k,r) = V>(k)A(r) 



(1) 



transformed according to either identity, (^(k)) 7^ 0, or 
non-identity, (i/)(k)) = 0, one-dimensional representation 
of the point group symmetry of the crystal—. Here and 
after, the angular brackets mean the averaging over the 
Fermi surface, ipfe) are the functions of irreducible rep- 
resentations, (|-0(k)| 2 ) = I. The generalization for the 
multidimensional superconducting states can be easily 
considered. 

The derivation for the case of superconducting state 
with an order parameter transforming according to gen- 
eral identity representation leads to quite cumbersome 
expression for the free energy. We shall consider the sim- 
plest example of identity representation with ip(k) = I, 



or s-wave pairing superconductivity, where the free en- 
ergy functional is 



F = J d 3 r l^a\A\ 2 + irN A* 
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Here, 



1 i/iH\ / 1 i[iH{) 

2 ~ 2^T ) ~ V 2 ~ 2ttT 
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(3) 



and Hq = Hq(T) is the critical field in homogeneous su- 
perconductor determined by the equation 



ln^=K 



The coefficients 



1 i/i-ffo 

2 ~ 2nT 



K nm = 2mJ2 
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'f=o (uj u - inH) n (Cj y - ifiH) m 
and oj v = irT(2v + I) are Matsubara frequencies, 

sign uj v 

LO v = LO v H , 

It 

D = — iV + (27t/$ )A, v(k) is the Fermi velocity, N is 
the density of states at the Fermi level. We put through 
the article Ti = 1. 

Near the critical temperature T c , both fields Hq(T) 
and the upper critical field tend to zero and the latter 
should be determined from the linear Ginzburg-Landau 
equation giving well known Gor'kov resulli^i with small 
correction due to paramagnetic effect. On the contrary, 
near the tricritical point, at large enough Maki parameter 
the upper critical field is close to Ho (T) such that one can 
write 



a = a (H - H ) 



2irT 



2-kT 



(0) 



and put H = Hq in all other terms of the functional. 

The free energy functional for non-identity representa- 
tion is 



F= I d 3 r{&\A\ 2 
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Here. 
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and i?o = Hq(T) is the critical field in homogeneous su- 
perconductor determined by the equation 



ln^=5R 




(9) 



Near the critical temperature T c , both fields H (T) and 
the upper critical field tend to zero and the latter should 
be determined from the linear Ginzburg-Landau equa- 
tion. Near the tricritical point, at large enough Maki 
parameter, the upper critical field is close to Hq(T) such 
that one can write 



a (H - H ) 

N /i(H - Hp] 
2ttT 



1 _ ijHIo 
\ 2 2ttT 



1 



4tttT 



(10) 



and put H = Hq in all other terms of the functional. 

In addition to the superconducting energy © or (7J), 
one should in principle also include the magnetic energy 



F„ 



1 

8^ 



d 3 r(H - H cxt ) 2 , 



(11) 



where H oxt is the external field, in the total free energy 
functional. In the following, we will negelect the con- 
tribution of such term by assuming that the screening 
current are not important (high-K limit) and H = H cxt . 



III. CRITERIA FOR THE APPEARENCE OF 
FFLO STATE AND THE FIRST ORDER 
TRANSITION 

A. Identity representation 

With the purpose to derive simple analytic criteria for 
the appearence of FFLO state and the change of the sec- 
ond order normal metal - superconductor transition to 
the first order one, let us make the angular averaging in 
the expression for the free energy in the case of pure 
s-pairing (^(k) = 1) in a metal with spherical Fermi su- 
face: 
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where v is the modulus of Fermi velocity and the term 
A~ 4 = (2eH/c) 2 originates from noncommutativity of the 
operators D x and D y . This value serves as the measure 
of the orbital effect such that the orbital effects free sit- 
uation corresponds to the limit A — > 00. 

Let us choose the magnetic field direction along the z- 
axis A = (0, Hx, 0)). So, for the Abrikosov lattice ground 
state A = fo{x,y)f{z) which is the linear combination 
of Landau wave fuctions with n = multiplied by expo- 
nentially / exp = exp(iqz) or sinusoidally / S j n = \/2smqz 
modulated function along z-direction, one can substitute 
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Making use of the properties: 



(D ± A) 2 = 0, J rf 3 r|A| 2 |D x A| 2 = ± J d 3 r|A| 4 , 



we come to the free energy in the following form 



(14) 



F = J dxdy {a (H - H{q)) |<^ | 2 + ixN B\ Vo \ A } , 
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H(q) = H a 
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and the coefficient B is given by 



B = 



K 30 
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in the conventional superconducting vortex-lattice state, 

^ = *-»V(^-t?) (18) 

in the exponentially modulated FFLO phase, and 

Wo q 2 v 2 (MUi K i2 \ 
B sm = - 2 r ^_-_j (19) 

in the sinusoidally modulated FFLO phase. 

The critical field value H c is found by taking the max- 
imum of H = H(q) as the function of q. The usual 
superconducting state appears at q = 0, while the FFLO 
state is formed when the maximum of H is reached at 
finite q = qo, where 

10 T 36t 

The FFLO state appears when the coefficient at q 2 in the 
square brackets of Eq. I|16|) changes the sign from positive 
to negative and it exists at 

Let us now examine the question of the transition type. 
It is determined by the sign of the coefficient at \<fio\ 4 in 
the expression (jlBfl for the free energy. Hence, the first 
order transition occurs at B < for the transition into 
usual superconducting state. It occurs at -B e xp(<7o) < 
for transition into FFLO state with exponential modula- 
tion or -B s i n (go) < for transition into FFLO state with 
sinusoidal modulation. 

To see explicitly the role of orbital effects and the im- 
purities in the formation of FFLO state and the change 
of the transition type, let us look on them separetely. 

In the clean, paramagnetic limit (A, r = oo), Eq. 1151) 
coincides with the free energy derived in Ref^. There, 
the inequality -K30 < was obtained as the condition 
both for the change of the transition type from normal 
to uniform superconducting state and for the FFLO state 
formation. In the vicinity of the tricritical point, keep- 
ing H = Hq(T) in ivjo, one finds that K30 changes its 
sign as a function of the temperature when T = T* . 
Therefore, K^o oc (T — T*) has the meaning of an ef- 
fective temperature close to this point. Further study 
reveals that, while B c ^ p (qo) — — .K30/6 remains positive 
at T < T t , B sm (qo) = A' 30 /36 becomes negative. This 
means that the first order transition from sinusoidally 
modulated FFLO state is favored at T < T*. The quali- 
tative superconducting phase diagram which results from 
this study is shown on Fig.QJa). 




temperature 



FIG. 1: Qualitative superconducting phase diagram for 
three-dimensional s-wave superconductor in presence of 
strong paramagnetic effect. The clean, purely paramagnetic 
case is shown on (a) , the purely paramagnetic case with disor- 
der is shown on (b), and the clean case, with some orbital ef- 
fect is shown on (c). Possible phases are the normal state (N), 
the conventional superconducting state: uniform state (U) in 
the absence or orbital effect, Abrikosov vortex lattice state 
(A) in presence of orbital effect, and the FFLO modulated 
state with exponential modulation (FF) or sinusoidal modu- 
lation (LO) along the applied field, eventually with Abrikosov 
vortex lattice (A+FF or A+LO) if orbital effect is present. 
Thin lines correspond to second order transitions, thick lines 
correspond to first order transitions, dashed lines correspond 
to transitions between different superconducting states and 
they have not been calculated in the present work. 
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In presence of impurities but neglecting the orbital 
effect (A — * oo), we obtain from Eq. (|21[1 inequality 
K21 < as condition of FFLO formation. It is easy 
to check that at the temperature determined by equa- 
tion K21 = 0, where the coefficient K21 changes the sign 
and, hence, the finite q modulation appears, the coeffi- 
cient K30 is already negative. Therefore, Bq is negative 
and the normal metal transforms to homogeneous super- 
conducting state by means of the first order transition. 
The impurities shift FFLO state to lower temperatures 
leaving unchanged the temperature of change of type of 
transition. The qualitative phase diagram in this limit is 
shown on Fig.QIb). 

On the other hand, taking into account only orbital 
effects, that is in the completely pure case (r — ► 00), one 
can rewrite the condition l|21[) of FFLO appcarence as 



K30 < 



3 v 2 K, 



50 



10 A 2 



(22) 



In pure case and at A — > 00, the FFLO state appears ex- 
actly when the coefficient -K30 changes its sign. Whereas 
in presence of orbital effect, the FFLO state appears at 
slightly negative K30 determined by negative value of 
K$q. Moreover, the condition of the first order transition 
into usual superconducting state, Bq < 0, is rewritten as 



A" 



:SII 



I v 2 K 5Q 
3 A 2 



(23) 



The comparison of these two inequalities makes clear 
that, due to the orbital effect, the change of the type 
of transition always appears at lower temperature than 
the FFLO state formation. This conclusion is in cor- 
respondence with the results of the pape»A& where the 
qualitative phase diagram shown on Fig. ^c) was first 
proposed. 

Thus, we find that impurities and orbital effect act in 
opposite directions regarding the shift of the tempera- 
ture below which FFLO state will appear. In the fol- 
lowing, we study the interplay between low impurities 
and orbital effect (1/A, 1/t — > 0) on FFLO state for- 
mation. In this limit, the change of the transition type 
at the normal/conventional superconducting vortex lat- 
tice transition, determined by Bq < 0, is still given by 
Eq. (|23() . The temperature below which transition into 
FFLO state occurs is determined by Eq. I|2I(I . In leading 
order in I/A 2 and 1/r, this equation yields: 

1 K A0 3 v 2 K 50 
K ™<2~ + I5^^> ( 24 ) 

where we made use of the property : 

m 
It 

In Eqs. (ESJl, HI, we recall that K 30 oc (T — T*) 
has the meaning of an effective temperature, while K^q 
and K 50 have to be evaluated at the tricritical point 
(T* , H*), where they take negative values. By comparing 



Eqs. (JJlJl, we find that the FFLO state appears at 
temperatures higher than the critical temperature when 
the impurity concentration remains low enough: 



K. 
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t < I5A 2 K. 
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(25) 



When Eq. I|25|l is obeyed, the free energy l|15fl also al- 
lows to discuss the structure of the FFLO state which 
is realized at the second order normal/FFLO transition. 
When -Bexp(go) < B s - ln {qQ) (and both of them positive), 
with go given by Eq. 12U|) . the FFLO state with expo- 
nential modulation is energetically favored. In the limit 
of low impurity and orbital effect that we consider, this 
inequality corresponds to: 



K30 > 
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The FFLO state with sinusoidal modulation is favored 
when -B s in(<?o) < -Bexp(q'o)- The transition into this state 



becomes of the first order when -B s in(<Zo) < 0j that is: 



K 30 < 



25 #40 3 v 2 K 50 
4 t 4 A 2 ■ 



(27) 



The above discussion is summarized on Fig. It rep- 
resents the nature and type of the transition from normal 
to the different superconducting states as temperature is 
lowered, for a given ratio between orbital and impurity 
effects. 

When the transition is of the second order, minimiza- 
tion of the free energy H15f) on ipQ yields that the vortex 
lattice is triangular^ 

In the Refi^i, it was predicted that a very large Maki 
parameter may favor not only FFLO modulation, but 
also an order parameter formed of higher-level Lan- 
dau functions: A ~ (p n (x,y)f(z), with n > 0. In 
Eqs. . lfH)|) the expression (A~ 2 + q 2 ) should be then 
replaced by ((2n + 1) / X 2 + q 2 ) , and the fourth order term 
B in the free energy l|13|) should also be calculated ac- 
cordingly. In particular, the condition to maximize the 
critical field l|16fl at the second order transition from nor- 
mal to vortex lattice state with Landau functions of level 
n > 1 is: 



K <r 1 Km 
< 2 — 



3(2n + 1) v 2 I<! 



10 
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(28) 



This equation is obtained by requiring that g 2 > in 
Eq. J20J after substituting 1/A 2 by (2n+l)/A 2 . One can 
note that, when such inequality is obeyed, the transition 
from normal to vortex lattice state (with n = 0) has al- 
ready turned from second to first order, both in the "low 
impurity" case in the presence of sinusoidal modulation 
(because Eq. I|27|) is already obeyed), and in the "high 
impurity" case (because Eq. l(2^ is already obeyed). 

Therefore, we get now the qualitative picture of the 
superconducting phase diagram in three-dimensional s- 
wave superconductors with strong paramagnetic effect. 
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FIG. 2: This figure shows the nature of the supercon- 
ducting state which is realized just below the upper crit- 
ical line, as well as the type of the transition, for three- 
dimensional s-wave superconductors with spherical Fermi sur- 
face, at large paramagnetic, small orbital and impurity ef- 
fects, and for temperatures close to the tricritical temperature 
T* . X 2 K 30 /v 2 \K 50 \ cc (T - T*) is the effective temperature, 
(xt is the ratio between orbital and im- 
purity effects. N/A indicates the place where the transition is 
from normal state into conventional Abrikosov superconduct- 
ing vortex lattice state, N/A+FF indicates the place where 
the transition is from normal state into FFLO state with vor- 
tex structure in the plane parpendicular to the magnetic field 
and exponential modulation along the field, N/A+LO indi- 
cates the place where the transition is into FFLO state with 
vortex structure in the plane perpendicular to the magnetic 
field and sinusoidal modulation along the field; in the white 
(gray) region, the transition into superconducting state is of 
the second (first) order. 



At large impurity concentration, the transition from nor- 
mal to usual superconducting state becomes of the first 
order at low temperatures, while FFLO state may ex- 
ist at even more lower temperatures either as stable or 
as metastable state. On the other hand, at low im- 
purity concentration, while temperature is lowered, the 
phase diagram shows the second order transition from 
normal to usual supercondutcing state, then to exponen- 
tial FFLO state, then to sinusoidal FFLO state, and fi- 
nally the change of the transition order into such state. 
These conclusions are summarized in the phase diagrams 
shown in Fig.n an d[21 



B. Non-identity representation 

As an example of similar calculations for nonconven- 
tional superconductivity we consider the d-wave super- 
conducting state ip(k) oc k 2 — k 2 in tetragonal crystal 
under magnetic field along c-axis (z-direction). One can 
rewrite first the Eq. Q) in the following form 



F= / d 3 r{&\A\ 



itNqA* 



Kt 



I):', 



«M 2 ui/2)Di 



2 <)Dl) 



-^((H 2 (vD ± ) 4 )+3(HM- 2 2 >Di^ 2 



MDt) 



■g|(^i/2)Di + (^)^| 2 



+ttN 



03 



#04 

8r 



IAI- 



4 «I/2)(DxA)*D x A 



-A^|A| 2 
+ (\4^v 2 z ){D z A)*D z A 
+7rA r x^| A ,2 



16r 



! W i/2)(DxA)*DxA 



|2„,2 



v 2 z )(D z A)*D z A 



+ (. 

+^V ff [(A*) 2 ((HM/2)(DxA) 2 
H^f{Wvl){D z Af)+c.c] 
-*M>§* [(A*) 2 «HM/2}(DxA) 2 
+(A*) 2 (H\ 2 v 2 z )(D z A) 2 ) + c.c.}}. 



(29) 



From this step, unlike to the conventional supercon- 
ductivity, the continuation of calculation for d-pairing in 
closed analytical form is not possible. The point is that 



the average (vDx) 



l (D+v- +D-v+) 4 ) con- 



tains the terms (|V>| 2 (-D ± i' =F ) 4 ) which are not equal to 
zero in tetragonal crystal. Here D ± = (D x ± iD y )/\/2 1 
V = (v x ± iv y )/\/2 and v\ = v 2 + v 2 . Hence, un- 
like to conventional superconductivity, the Abrikosov lat- 
tice ground state in tetragonal superconductor with d- 
pairing^ is the linear combination of functions consist- 
ing of infinite series of Landau wave functions (p n (x,y) 
with n=0, 4, 8, 12. . . multiplied by exponentially or sinu- 
soidally modulated function f(z) along z-direction 



A = f(z)(A tp + + A$tp 8 ....) 



(30) 



Fortunately, in the limit of large Maki parameters we arc 
interested in, one can work with cut-off series of the form 
similar to s-wave pairing 



A « f{z)<p 



(31) 



and also neglect the terms like (\ip\ 2 (D^ v^) 4 ) in the 
Hamiltonian (the proof of this property is found in Ap- 
pendix C). So, we shall use the substitution 

<H 2 (vDx) 4 )A 

=> (\4>\ 2 (irv + ) 2 ){D-D + D-D + + D'D'D+D^A 
= 3(\ij\ 2 (v-v + j 2 )\- 4 A (32) 

Thus, the further calculations have the sence of varia- 
tional treatment. 
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Similar to the case of conventional superconductivity, 
we now obtain from Eq. (|29[l 



F = 
where 



dxdy {&o (h - H(q)) |^ | 2 + ^(O)^! 4 } , 

(33) 



H(q) = H - 



a 



#15 I / , 2 \ 1 
\\WV i ) =■ 

32r V ±, 2X 2 



and 



5 K03 .. . 4 . Kq4 #05,,, ,4 2 , 
Bo = — <H >__-^.(H „ ±< 



64rA ; 



(34) 



2„,2 



(35) 

for the conventional superconducting vortex lattice state, 



Bcxp(q) = B 



3#o 



■<HM X 



3A" 



IMi 



16r 



for the exponentially modulated FFLO state, and 



(36) 



B sin (q) 



3B , g 2 



or, 



4„,2\ 



7# 



0(> 



16r 



(MM 



(37) 



for the sinusoidally modulated FFLO state. 

The study is now similar to the previous section. 

The critical field H c is determined by the maximum of 
H{q) as the function of q. The FFLO state arises when 
the maximum of H occurs at finite wave vector 



2 
% 



(38) 



2„,2\ 3if 05 (|y>lM^) K la ft((4,'v 2 ± )(4wD) 



4A 2 



8rA 2 



2 4r 



The FFLO state appears with the sign change of the 
coefficient at q 2 in Eq. I)34|) and it exists at 



#03 < 



AX 2 (M 2 v 2 ) 



(39) 

The type of transition changes from the second to the 
first order with the sign change of the coefficient B at 
|(/3 | 4 in Eq. So, the first order transition from 

normal to Abrikosov vortex lattice state persists when 
Bq < 0, that is in the region of validity of the following 
inequality 



#03 < 



#, 



04 



5#06 (M 2 v 



M> 



2r 



2A 2 



16rA 2 (|V>| 4 ) 



(40) 



To see explicitly the role of orbital effects and the im- 
purities in the formation of FFLO state and the change 
of the transition type, let us look on them separetely. 
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clean, paramagnetic limit 



LO 



paramagnetic limit + disorder 




temperature 



FIG. 3: Qualitative superconducting phase diagram for quasi 
two-dimensional d-wave superconductor in presence of strong 
paramagnetic effect, (see also legend of Fig. 0. 



In pure limit r — > oo, the two inequalities l|39|l and 
(I40|l take the much simpler form, where FFLO exists at 



#03 < 



3tf„6<MVL«2> 



4A 2 (|^| 2 u 2 ) 

and the first order transition to vortex lattice state at 



(41) 



#03 < 



(mm: 

2A 2 (M 4 } 



(42) 



In quasi two-dimensional case, we deal in first approxi- 
mation with cylindrical Fermi surface: ipk = \/2(k 2 — k 2 ) 
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and v± ~ v is constant. Then, the inequalities look even 
simpler: FFLO state exists at 



K 03 < 



3K 05 v 2 
4A 2 



(43) 



and the first order transition to vortex lattice state at 

K 05 v 2 



K Q3 < 



2A 2 



(44) 



The comparison of these two expressions makes evident 
that the critical point where the second order transition 
from normal metal to superconductor transforms to the 
first order one lies at higher temperature than the point 
at which FFLO state arises. The resulting phase diagram 
is qualitatively shown on Fig.|3Jc). It corresponds to the 
experimental observation in CeCoIn^^. 

This conclusion is just the opposite to the case of s- 
wave superconductivity in isotropic metal considered in 
previous subsection. It is obvious that this difference has 
pure geometrical origin (order parameter and the Fermi 
surface anisotropy) and does not originate from the dif- 
ference in applied theoretical approaches: gradient ex- 
pansion ini& and nonpertubative treatments^!. 

In absence of orbital effect (A — > oo), the FFLO state 
exists at K03 < 0. In the pure limit (r — > 00), it gives 
therefore higher critical field than the critical field cor- 
responding to the transition from normal to uniform su- 
perconducting state which also changes its type at the 
same place. In this pure, paramagnetic limit, K 03 = K30 
and Hq{T) = Hq(T). Therefore, the tricriical point de- 
fined by K30 — along the critical line H$(T) is still 
(T*,H*) for d-parining. The actual structure of the 
FFLO state which is realized just below the critical line, 
at T < T*, was considered in Refsi 17 i 18 i 20 i 26 for quasi 
two-dimensional d-wave superconductors. It was shown 
that, close to the tricritical point, the second order tran- 
sition takes place from normal state to sinusoidally mod- 
ulated state, with the direction of modulation parallel to 
the conducting planes and along the nodes of the order 
parameter. The resulting phase diagram is illustrated in 
Fig. Ela). In our analysis, we considered the modulation 
along the applied field. It yields B oxp (qo) = —3X03/8 and 
B B m(q ) = -K 03 /16, that is B exp (q ) > B Bm (q ) > in 
the region of existence of FFLO state. Therefore, we ob- 
tain the same topology for the phase diagram as shown 
in Fig. Ha). 

In absence of orbital effect (A — ► 00), but with some 
disorder, one can check that, at the point where FFLO 
instability takes place, defined by K03 = 0, Bq is still 
positive. Therefore, the phase transition from normal to 
uniform superconducting state is still of the second or- 
der. Thus the influence of impurities is also opposite to 
the case of s-pairing in isotropic metal. If we consider the 
modulation along the applied field, we find -B e xp(<7o) = 
-3X 03 /8 -A' 04 /8r and 4in(<Zo) = -K 03 /16- 3K 04 /16t. 
Thus, both are positive, due to negative value of -K04 at 
the tricritical point and transition from normal to FFLO 



state remains of the second type, while the sinusoidal 
(exponential) modulation is favored at K$3 < Kq 4 /5t 
(K 04 /5t < K 03 < 0). The corresponding phase diagram 
is qualitatively shown in Fig. E{b) . It has the same topol- 
ogy as the phase diagram given for modulation along the 
nodes of the order parameter in the presence of impurities 
in Ref,2&. 

When both orbital effect and impurities arc present 
and small, and assuming an almost cylindrical Fermi sur- 
face, we can obtain the more precise following picture for 
the phase diagram. In the leading order in 1/A 2 , 1/r, one 
finds that the equations defined by l|39|) and l|4U|l cross at 
K a4 /2r = 3v 2 K a5 /8X 2 . When 



I 3v 2 K 05 
t 4A 2 A' 04 ; 



(45) 



transition from normal to usual superconducting vortex 
state turns from second to first order when tempera- 
ture is lowered. The modulated FFLO state may ex- 
ist as a stable or metastable state below the first or- 
der transition line, with eventually even larger critical 
field at lower temperatures (Fig. EJc)). On the other 
hand, when inequality (|45|l is reversed, the impurity ef- 
fect dominates: as the temperature is lowered, the tran- 
sition from normal to superconducting state is changed 
into the transition from normal to exponentially modu- 
lated state at K 03 < 3v 2 Ko 5 /4X 2 . Then, it is changed 
into the transition into sinusoidally modulated state at 
Boxp(qo) > B sin (q ). With Eqs. JHJ), (23), JSJ, and 
appropriate averaging over Fermi surface, one finds that 
this occurs at: 



K03 < 



K 



04 

5r 



3v 2 K 



05 



5A 2 



(46) 



In this region, B cxp (qo) and B s i n (qo) remain positive, 
therefore the critical line remains of the second order. 
These results are summarized in the Fig. QJ 

From the previous discussion, one can guess that our 
Ansatz <|33[1 is not the most general. Indeed, in the pres- 
ence of orbital effect, one should also consider order pa- 
rameters in the form 

A = f(z)(A n if n + A n+4 Lp n+4 + ...), 

where n > 0. Therefore, the FFLO phases illustrated in 
Figs. 13 and 01 may compete with phases corresponding to 
order parameter with higher Landau level n > 0. This 
problem is reserved for further study. 

Note that these results are very sensitive to the shape 
of the Fermi surface. In particular, different topology 
for the phase diagram would be obtained for anisotropic 
d-wavc superconductor with elliptic Fermi surface. 



IV. FFLO INSTABILITY IN DISORDERED 
S-WAVE SUPERCONDUCTOR 

According to Ref>i£, instability toward FFLO state for- 
mation is always present in s-wave superconductors in 
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temperature 



_ 4 




FIG. 4: This figure shows the nature of the superconducting 
state which is realized just below the upper critical line, as 
well as the type of the transition, for quasi two-dimensional 
d-wave superconductors, at large paramagnetic, small orbital 
and impurity effects, and for temperatures close to the tricrit- 
ical temperature T*. (See also the legend of Fig.|2J. 



the paramagnetic limit. In particular, assuming that 
the transition from normal to FFLO state is of the 
second type, such instability occurs at T* = 0.56T C 
in clean systems and at vanishingly small temperature 



Tr* 



-Ao/(21nrAo) in dirty ones. 



On the other hand, a large orbital effect is detrimental 
to FFLO instability, as shown in clean s-wave supercon- 
ductors in RcLM. Indeed, it was shown there that the 
FFLO instability only takes place when the paramagnetic 
effect, characterized by a Maki parameter a° M > 1.8, is 
strong enough. 

In this section, we adress the question of FFLO in- 
stability in s-wave superconductors with disorder. In 
particular, we find that the second order transition to- 
ward FFLO state exists for any disorder provided that 
the orbital effect remains small enough. In the dirty 
limit (rAo <C 1), the Maki parameter characterizing its 
strength must be very large: aj^ > — 1/(tAq InrAo). 

As we discussed in the previous sections, transition 
from normal to conventional superconducting state may 
change its type. In particular, in dirty systems, such 
change of the transition type was shown to take place 
below some critical temperature as soon as af f > 
l*im Therefore, the FFLO state may either exist as a 
metastablc state below the first order transition line into 
conventional superconducting state or take place by first 
order transition with even higher critical field. We do 
not discuss the question of the type of the transition in 
the present section. 

Let us now derive the result. 

At the second order transition in superconducting 



state, the linearized self-consistency equation IjBlll is 



A In — + 2ttT3? 

r 



7 (A- 



A = 0, 



where the differential operator {CJ) is: 



(A, 



da 



i 



4tt a, + ix±i 



(47) 



(48) 



The most general form of the solution for the gap at the 
second order transition is A = <po(x, y)e lqz , where q is the 
FFLO modulation vector and <p is the Abrikosov vortex 
lattice formed of lowest level Landau functions. Using the 
identity 1/X = / °° dse~ sX and the properties of Landau 
functions, we find that the operator {C v ) applied to A 
yields the eigenvalue: 



x exp 



ds exp(— stl v ) / du 
o Jo 

-s 2 (l - u 2 ) cos f 



(49) 



8$n 



V 2 J 



At q = 0, Eqs. |I? )l .l|l9" | yield the second order crit- 
ical line H c2 (T) for transition into usual superconduct- 
ing vortex lattice state. In particular, in the paramag- 

(2) 

netic limit the upper critical field is H p = Ao/(2/x) = 
Hp/y/2 and it does not depend on the disorder. In the 
clean, orbital limit, the upper critical field is H^ 2 q = 
(7e 2 /27r)<f>oA§/u|,. At finite temperature and/or inter- 
mediate disorder, H C 2 must be found numerically 

In the dirty limit, the equations determining the upper 
critical field simplify greatly. Indeed, integration on s in 
Eq. I|49|) is cut off by impurity scattering time, s < r. 
Thus, assuming that H -C §o/(vft) 2 (as can be checked 
consistently later), we may expand the second exponen- 
tial in Eq. I|49|l and perform the integration explicitly: 



(50) 



Thus, we obtain the implicit equation for H C 2 in the dirty 
limit2i: 



T 

1b 5 t = » 



- * 



1 . .^H c2 
h l 

2 2ttT 



DH, 



<:2 



2$nT 



(51) 

where D = Vpr/3 is the diffusion constant. In particular, 
the critical field at zero temperature, 



H c2 (0) 



(2) 

interpolates between Hp in the paramagnetic limit and 
the upper critical field in the orbital, dirty limit, H^ 2Q = 
<S>oT c /2^fD, with the Maki parameter in the dirty limit 
defined as: at, = y/2H^ 20 /H v = ^q/ttD. 



Tjd 
H c2Q 



1 



d2 



(52) 
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In general, the critical field defined by Eq. (|47|l also 
depends on g: H = H(q). The actual critical field cor- 
responds to the maximal value of H(q) with respect to 
q. When it is obtained for q ^ 0, second order tran- 
sition into FFLO state is realized. Along the critical 
line H(q = 0) at given impurity rate and Maki param- 
eter, the triple point below which such transition may 
occur is defined by dH/d(q 2 )\ q= o = 0. (One could check 
that dH/dq\ q= o — is always true.) In order to obtain 
dH/d(q 2 )\ q=0 = 0, one can expand Eq. (|T7|) up to the 
second order in q, and obtain: 



H 
A 



H c2 (T)+kq 2 - 
(2ttT) 2 



0(q 4 ) 



d(C v (q))/d(q 2 ) 



(53) 



q=0 



where ( = l/a d M + i. Making use of Eq. ij4*§|) and 
integration by part, one finally obtains the condition 
dH/d{q 2 )\ q=0 oc A = in the form: 



= SR 



(54) 



In particular, at zero temperature this equation defines 
the minimal Maki parameter which allows existence of 
FFLO state. In the dirty limit, Eq. (|5^jl is easily inte- 
grated at zero temperature: 



= 3^/ duj 



1 



i + K/) 2 ( i 

—T - t-Aq In 



a d M (rA ) 2 \a d M 



(55) 



It yields the critical Maki parameter above which FFLO 
instability exists: 



1 



rA In -4- 

u tAo 



(56) 



Therefore, FFLO instability is present in dirty supercon- 
ductors provided that the orbital effect is small enough. 



V. CONCLUSION 

In conclusion, we derived microscopically the gener- 
alized Ginzburg-Landau free energy functional which is 
adequate to describe conventional and unconventional 
singlet superconductors in the presence of paramag- 
netic, orbital and impurity effects. This free energy 
was used to predict the superconducting phase diagrams 
of three-dimensional s-wave superconductors and quasi- 
two-dimensional d-wave superconductors under magnetic 
field perpendicular to the conducting layers. These phase 
diagrams prove to be quite different and to be very sensi- 
tive to geometrical effects such as the nature of the order 



parameter and the shape of the Fermi surface. In par- 
ticular, we found that impurities tend to favor the tran- 
sition from normal state to the Abrikosov vortex lattice 
state, with the change of the transition type as the tem- 
perature is lowered in the s-wave case, while they tend 
to favor the transition from normal state to the FFLO 
state with vortex lattice plus additional modulation of 
the order parameter along the field direction in the d- 
wave case. We also found that the orbital effect acts in 
the opposite direction. That is, it tends to favor transi- 
tion from normal to the FFLO state in s-wave case, while 
it tends to favor the transition from normal state to the 
vortex lattice state with the change of the transition type 
in the d-wave case. 

In addition, we determined the criterion for instabil- 
ity toward non uniform superconducting state in s-wave 
superconductors in the dirty limit. 



APPENDIX A: FREE ENERGY FUNCTIONAL IN 
SUPERCONDUCTOR DOPED BY IMPURITIES 

The free energy functional expanded over the order 
parameter for superconducting state with pairing inter- 
action 

V(k,k') = -Vip{k)ip(k') (Al) 
has the following form 

z_ v 

q 

T E E G-(k,k')G--(-k + q,-k' + q')A£ !q A k ,. q , 
oj kk'q' 

+?£ E E G "( k >p) A U 

uj qi— q2=q4— q.3 kplm 

G- w (-k + qi ,-l + qi )A liq2 

G"(m - qi + q 2 , 1 - q x + q 2 ) A^ q3 

G-"(-m + q 4 ,-p + q 4 )Ap, q4 . (A2) 

Here, the order parameter is given by Fourier transfor- 
mation of Eqn 

A k: q = ^(k)A q = J d 3 rcx P Hqr)A(k,r), 

and G"(k, p) is exact electron Green function in normal 
metal with arbitrary configuration of impurities. Aver- 
aging of the free energy over impurity configurations de- 
mands calculation of averages of vertices^ 

^k,q = 

G "( k , k')G-"(-k + q, -k' + q')A k ,, q , (A3) 

k'q' 

obeying of equation 

A k , q = J2 G w (k,k') G-"(-k + q,-k'+q') 

k'q' 
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[Ak/, q / + nu 2 ^A p , q ,], 



(A4) 



where nu 2 = 1/2ttNqt, n is impurity concentration, it 
is the amplitude of scattering, and r is mean free time 
of scattering of quasiparticlcs. Substituting the Green 
functions by its average 



G"(k,p) = G w (k)<J(k-p), 
1 



G ±w (k) 



(A5) 
(A6) 



±<(w w - - f(k) 

we obtain from (|A4fl 

J 4 k:q = G"(k)G- w (-k + q)A k , q , (A7) 

where 

(A8) 



A k ,q = A k , q + JiU 2 ^ ^4p !q , 
P 



and 



EC"(k)G-"(-k + q)A k?l 



E A ^ = r 



- nu 2 G"(k)G-"(-k + q) ' 



(A9) 



Then, following the procedure developed in£i, after the 
averaging of free energy l|A2|l we obtain 



F = F 2 + F± 



(A10) 



where 



T££G^(k)G-"(-k + q)A£ >q A k , q ( A11 ) 

w kq 

^ 4 =|e e (e^wau 

w qi-q2=q4-q.3 I k 



qi-q2=q4-q: 
G""(-k + qi )A k , q2 G"(k qi + q 2 )A k 
G-"(-k + q 4 )A kiq4 

W^G-(k)A kjqi 

k 

G—(-k + qi )A u , q4 G-"(-k + q 4 ) 
]T G" w (-k + qi )A k , q2 G^(-k + q 4 ) 

k 

A kiq3 G-(k-q 1 +q 2 ) 

+nu 2 G'^i-k + q 4 ) A kiq4 

k 

G"(k)A k!q3 G w (k-q 1+ q 2 ) 
^G-(-k + qi )A k)qi 

k 

G"(k)A k , q2 G w (k- qi + q 2 )} 



(A12) 



The further calculations are different for different su- 
perconducting states. Expanding the quadratic in re- 
spect of the order parameter terms up to forth order in 
respect of q and performing £ integration we obtain for 
the case of s-pairing 

^G-(k)G-(-k + q)A k , q A k , q 



7rA n 



o U <(qv) 2 ) ((qvr) 



(((qy) 2 )) 2 

32x^3 



1 2 



(A13) 



Corresponding expression for superconducting order pa- 
rameter transforming according to non-identity represen- 
tation even in respect of k is 



£G"(k)G-(-k + q)A k)q A k , c 



ttA 



1 - 



((qv) 2 |^(k)| 2 ) , ((qv) 4 |^(k)P 



n„ { mi 
l((qv) 2 ^(k))| 2. 



i6m 



A q | 2 . 



(A14) 



Here 



il„ = uj v — iuH, il„ = u) v — ifiH, a>„ = u) v 



sign u> u 
(A15) 

Substituting (|A13f) and (|A14|I in ilAllfl and performing 
Fourier transformation to the coordinate space (accom- 
panied by substitution q — > D = — iV + (2-7r/<I>o)A) we 
come to the quadratic terms in the Eqns J5Jl and Q) cor- 
respondingly. 

In terms of the forth order in respect of the order pa- 
rameter one should calculate A kq up to the second order 
in q. For the s-pairing state it is 



and for a non-identity representation 



(A16) 



A 



i= U)_((sv)!t0U c 

' q \ V ; SrQ^ 2 i 



(A17) 



It is easy to check that in the latter case in terms of the 
forth order in respect of the order parameter and up to 
the second order in q one can put just A kq = ^(k)A q . 

Substituting l|A16|) and 1A17|) in (|A12jl . expanding the 
Green functions up to the second order in q, performing 
the integration over £ and Fourier transformation to the 
coordinate space we come to the quartic terms in the 
Eqns J2J) and Q correspondingly. 
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APPENDIX B: FREE ENERGY FUNCTIONAL 
DERIVED FROM EILENBERGER EQUATIONS 

In this appendix we propose another method to derive 
the superconducting free energies an d which are 
introduced in Sec. [H] 

The quasiclassical theory of superconductivity forms a 
convenient framework to study conventional and uncon- 
ventional superconductors in the presence of magnetic 
fields or impurities^ In this theory, the superconducting 
gap is related to the anomalous function /„(k, r) through: 



A(r)= 7 rTF5J^(k)*/,(k,r)) £) 



(Bl) 



where the brackets stand for the averaging over the Fermi 
surface labelled by k, V and V'(k) define the pairing in- 
teraction IjAljl . The anomalous function is determined 
by the set of Eilenbcrgcr equations: 

l -vT>f v + (£l v + ±-{g v ))f v = UA+^{f v )^g v , (B2) 

where 



g v = sign(w„)y 1 - f v ft and fl = -f*_ {v+1) . (B3) 

The magnetic field H = rotA is combined with gradient 
in D = — iV + (27r/$o)A, and f2„ — lu u — i/iH where io v 
is a Matsubara frequency. 

Near the second order transition from normal to super- 
conducting state, the order parameter A is vanishingly 
small. Moreover, we assume that the order parameter 
is slowly varying on the scale of the superconducting co- 
herence length. Then we may expand the selfconsistency 
equation i|Bl|) up to the third order terms in the gap, and 
the fourth order terms in the gradient expansion. 

In the following we proceed separately for the cases of 
identity and non-identity representation. 



where £„ = (f2„+ivD/2) Expanding up to the fourth 
order terms in the gradient expansion, we find: 

(1) = A ((vD) 2 )A ((vD) 4 )A ((vD) 2 ) 2 A 
w„ 4u>*£l v 16<D2f23 32™3f]3 ' 

(B6) 

where oj v = cj„ + l/2r. 

In the second order in A, we find from Eq. I)B3|1 : g„ = 

In the third order in A, on gets from Eq. I|B2|I that 

(3) 

fi> is the solution of the linear differential equation: 



lvD/( 3 )+O i/ /( 3 )=A 5 ( 2 ) + i(/0 



(3)\ 



2t 



(B7) 



(3) 

We can now express (f v ) in terms of A and we make an 
expansion up to terms of the second order in the gradient 
expansion. We obtain: 



(3)\ 



AIAI 



1 



[<(vD) 2 )A|A| 



2^ scdini 

+ ((vD)(|A| 2 (vD)A + A 2 (vD)*A*)) 
-A(|vDA)| 2 ) + |A| 2 ((vD) 2 A) + A 2 ((vD*) 2 A*)] 

[((vD) 2 )A|A| 2 + ((vD)A 2 (vD)*A*) 



16™^2 



+ |A| 2 ((vD) 2 )A + A 2 ((vD*) 2 )A*)] . 

(B8) 

At ui v of arbitrary sign, one should substitute Q v — > 
\uj v \ — i/ii?sigiu<j„ and £7„ — > \u v \ + l/2r — i/i_Hsigncj„ in 
Eqs. CBS), Est . 

Inserting now Eqs. (|B6() . (|B8() into (|B1|I . we can put 
the selfconsistency equation for the gap in the form: 



1. Identity representation 

For simplicity, we only consider identity rcprescnation 
with ip = 1. 

The set of Eilenberger equations can be expanded per- 
turbatively in the gap. In this expansion, g v only con- 
tains even terms: g v = gi°^ + g v 2 ^ + . . ., while /„ only 
contains odd terms: f v = f„ + f^ + . . .. 

In the zeroth order in A, at u u > 0, we find gl 0) = 1. 

In the first order in A, is the solution of the lin- 



earized differential equation (|B2(I : 

l v D/«+^/«=A + ^(/«), (B4) 

where VL V = Q u + l/2r. Expressing (fv ) in terms of A, 
we get: 



(fW) = A 



(B5) 



= ( m _L + 27rT £ I _ nKw A + ^i((vD) 2 ) A 
-^p ( (vD)> - ^((vD) 2 ) 2 A + ^A|A| 2 



16r 



[4|A| 2 ((vD) 2 )A + A 2 ((vD*) 2 )A* 
+3A*((vDA)) 2 -2A|((vDA))| 2 ] 



3|A| 2 ((vD) 2 )A + A 2 ((vD*) 2 )A* 
+2A*((vDA)) 2 - 2A|(vDA)| : 



(B9) 



where we used the standard regularization rule: 



UJ>0 



and the coefficients K nm are defined in Eq. JSJ). We can 
check straightforwardly that Eq. (|B9|) corresponds to the 
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saddle point equation for the free energy functional J2J: 
SF 



5A*(r) 



= 0. 



(BH) 



2. Non-identity representation 

One should proceed along the same line to derive the 
gap equation for non-identity representation (when (ip) = 
0). In particular, one can obtain: 



(« iJ ) = ~- 



A (|V>| 2 (vD) 2 )A (|Vf(vD) 4 )A 



n„ 4^3 16^5 

| (f(vD) 2 )(y;(vD) 2 )A 



(B12) 



and 



* f(3)\ 



i 

2&, 



AIAI 



+ [(vD) 2 A|A| 2 + vD (|A| 2 (vDA) 

+A 2 (vDA)*) + |A| 2 (vD) 2 A 
+A 2 (vD 2 A)* + A|vDA| 2 ]) 

-^-(|^| 2 [(vD)| A| 2 (vD)A + (vD) 2 |A| 2 A 
16riZ° 

f 2|A| 2 (vD) 2 A + A 2 (vD*) 2 A* + A|vDA| 2 ]). 



(B13) 

Inserting now Eqs. l|BT2llBT3)l into JHTJl, we can put 

the sclfconsistcncy equation for the gap in the form: 







(ln-^ + 27rTVi-7rA- i) A 



*(vD) 2 )A 



16 



■(M 4 (vD) 4 }A 



[4|A| 2 ((vD) 2 )A + A 2 ((vD*) 2 )A* 
+3A*((vDA)) 2 - 2A|((vDA))| 2 l 



APPENDIX C: LOWEST LANDAU LEVEL 
APPROXIMATION 

In this Appendix we shall prove that in the limit of 
small influence of orbital effect (large Maki parameters) 
the function given by Eq. (|31|l is appropriate variational 
function for the Abrikosov lattice ground state in tetrag- 
onal superconductor under magnetic field directed along 
c-axis. With this purpose let us consider the Hamiltonian 
of the form 



where 



H = H Q + Hi, (CI) 



Ho = aa~a + f3[a~a a~a + + a~ a~ ' a+a 
+a~ a + a + a~ + a + a~a~a + 

+a + a + a~a~ + a + a~a + a~], (C2) 



and 



Hi = 7 (a + r + 5{a-y 



(C3) 



The dimensionless differential operators a = AD act 
on the Landau states 4> n (x, y), n = 0, 1, 2, ... as follows 



a 4> n = Vnfin-i, a + 4> n = y/n+ l(f> n+ i, (C4) 
such that 

Ho&o = £o0o, £o — a + 3/3, (C5) 

and 

H cj) A = e 4 04, £4 = 5a + 123/3. (C6) 
Let us consider variational wave function 

4>= 4>Q+ O404 (C7) 

with a4 as a variational parameter and calculate the ex- 
pectation value 



/ dxdy<pH(f) _ e + V2A(j + S)a 4 + 
J dxdy\4>\ 2 



(C8) 



The minimum of this expression is determined as solution 
of the equation 



2(e 4 - e ) 



a 4 -( 1 + S)=0 (C9) 



+ !* r5|A| 2 ((vD) 2 )A + 2A 2 ((vD*) 2 )A* 

16r 

+3A*((vDA)) 2 - 4A|(vDA)| 2 ] . 

(B14) 



It is clear that at 7 = 6 = in other words at H = Hq 
the variational parameter a 4 = 0. In general 



a A = v- v\j\ + l/v 2 as - — , 

2v 



We can check straightforwardly that Eq. (|B14|) corre- 
sponds to the saddle point equation for the free energy 
functionnal Q. 



where 



£4 — £q 



24(7 + 5) 



(CIO) 



(Cll) 
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The values of coefficients we used are 

nN K 03 (\^vl) 
a 8A 2 



Hence, v cx A 2 and our variational parameter proves to 
be small as 



(C12) 



and 



7 = S = 



64A 4 



7rN K 05 (\^(vi-6vlvl+v 4 y )) 
64A 4 



(C13) 



(C14) 



( |f 



where, for brevity we have written a, /?, 7 and 5 in clean 
limit r — > 00). Thus, in the limit of large A we obtain 



16A- 03 A 2 (HX)| 



VG(-K 05 )(\^(vi-6vX + vt))\ 



(C15) 



— . (C16) 
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